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Classical dynamics of a two-species 
Bose-Einstein condensate in the presence of 
nonlinear maser processes 



B. M. Rodriguez-Lara and R.-K. Lee 



Abstract The stability analysis of a generalized Dicke model, in the semi-classical 
limit, describing the interaction of a two-species Bose-Einstein condensate driven 
by a quantized field in the presence of Kerr and spontaneous parametric processes is 
presented. The transitions from Rabi to Josephson dynamics are identified depend- 
ing on the relative value of the involved parameters. Symmetry-breaking dynamics 
are shown for both types of coherent oscillations due to the quantized field and 
nonlinear optical processes. 



Qi 1 Introduction 



A great deal of attention has been granted to research about unifying concepts 
^ in classical and quantum physics through experimental demonstrations involving 

Bose-Einstein Condensates (BECs) at a macroscopic scale llj. For example, by 
considering BECs loaded in a double-well potential, the quantum tunnelling be- 
V/-^ tween two trapped condensates provides a possibility to study and understand 

^J symmetry-breaking, self-trapping, and Josephson oscillation; all of them are fun- 

(^ damental problems in quantum physics. Macroscopic quantum self-trapping |2 1 and 

psj bosonic Josephson junction |3| have been demonstrated recently, and successfully 

^-H described by a mean-field approach. The realization of these macroscopic quantum 

j>. self-trapped modes opens up new avenues for research; e.g., generation of squeezed 

• ^ atomic states L4J and atomic interferometry |i5j. 
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By using the symmetric and asymmetric stationary eigenstates of the Gross- 
Pitaevskii equation for a macroscopic condensate trapped in a symmetric double- 
well potential, the nonlinear Hamiltonian (in units of h) describing coherent atomic 
tunnelling between two zero-temperature BECs is equivalent to that of a nonrigid 
pendulum OIJl, i.e., 



A 



HsDW = '^^- ^y-z^ cos 0, 



(1) 



where the length of pendulum decreases with the angular momentum z, which de- 
notes the fractional relative population between the condensates in the two wells. 
The ratio of the on-site interaction energy and the coupling matrix element is char- 
acterized by the parameter A, while the tilt angle shows the phase difference 
between the two condensates. 

This nonlinear system allows for macroscopic quantum self-trapping {i.e., lo- 
calized oscillations of the fractional population difference z), which can be easily 
deduced from the equations of motion of the system. 



z = — V 1 — z^ sm<P 

z 



<P = Az- 



VT^- 



COS0. 



(2) 
(3) 



The fixed points of such equations of motion present a pitchfork bifurcation at the 
critical point A^ = 1 ; i.e., {zji = 0, = tt} for A < 1 and {zn = ±VA^— XjA , ^ = 
Tt} for A > 1. This bifurcation signals a symmetry-breaking in the dynamics of the 




Fig. 1 (Color online) a) Fixed points for the equations of motion of a BEC in a double-well 
potential, described by Eqs.ppl, at ^ = TT. The bifurcation of a single stable fixed point in the 
Rabi regime (solid purple) split into two fixed points in the Josephson regime (dotted blue and 
dashed red) is shown at the critical value, A^ = 1 ■ The unstable fixed point, or separatrix, in the 
Josephson regime (dash-dotted black) is also shown, b) Examples of trajectories for identical Rabi 
oscillations (solid purple) around the fixed point A (A = 0.5) and localized Josephson oscillations 
around the fixed points B (dotted blue ) and C (dashed red ) (A = 1.5) are demonstrated with 
symmetrical initial conditions {z± = ±0.5, 4> = Tt\. 
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system from Rabi to Josephson dynamics; such that for symmetric initial conditions, 
{z± — ztzQ,'P — Tt}, the trajectories will be identical in the Rabi regime, A < 1, but 
localize in opposite hemispheres of phase space in the Josephson regime, A > I. 
Figure [T] shows an example of such symmetry-breaking phenomena. The transi- 
tion from Rabi to Josephson dynamics due to the pitchfork bifurcation of the fixed 
points has been confirmed experimentally in an equivalent model consisting of two 
hyperfine states of a single atomic specie BEC coupled by a classical two-photon 
transition in the semi-classical limit ^. 

Recently, strong coupling between a BEC and the quantized field mode of an 
ultrahigh-finesse optical cavity has been demonstrated |9, 10, 11 1. The spectra of 
these strongly-coupled systems shows a level splitting attributed to different hyper- 
fine structures of the given atomic species. Motivated by these experiments, here, 
we consider a gaseous BEC, composed of bosonic atoms populating two hyperfine 
levels, coupled to a quantized field cavity mode, including nonlinear interactions 
among the ultracold atoms. Instead of two trapped BECs interacting through quan- 
tum tunnelling in a double-well potential, the hyperfine levels of the two-species 
condensate are driven by the quantized cavity field. It will be demonstrated later 
that the Hamiltonian describing this system consisting of a BEC in a cavity can be 
reduced into a nonlinear Dicke model with an additional atom-atom quadratic inter- 
action term; in the semi-classical limit, this system can be though of as a pendulum 
with changeable mass. In order to study a general model, we consider second and 
third order nonlinear processes for the field. Our goal is to present a steady-state 
analysis of the equations of motion in the large-ensemble-size limit for the collec- 
tive dynamics of this system. In Section l2j the model and possible physical real- 
izations are presented and discussed. In order to present analytical results, a weak 
regime is defined in Section l3] for weak coupling and nonlinearities (compared to 
the frequency of the driving field). The symmetry of Josephson dynamics is shown 
to break by the driving quantum field. Finally, Section|4]closes with a summary. 



2 Model 

The proposed model Hamiltonian describing the interaction of a single electromag- 
netic cavity mode a coupled to a BEC with two internal hyperfine structure levels /7| 
and /?! can be obtained from the Gross-Pitaevskii equation describing a two-species 
BEC interacting with a quantum field 1 12 1 and is written as. 



Ho = cooa^a+ £ {Efb]bj + \Gjfb]b]bfbj) 

{a + flt) {t\b^ + 1^^ + G^ib\blbJ,^, (4) 



+ -/A. 



where the frequency (Oq is the cavity mode frequency, the parameter g is the cou- 
pling strength between the cavity mode and condensed atoms, and A^ is the number 
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of atoms in the condensate. The energies of two internal hyperfine levels are la- 
belled as £1 and E-f with an intra-atomic transition frequency (Oa =E^ — E^. Here, 
we suppose the interaction Hamiltonian as a coupled two-component BEC fT3l, by 
introducing Gjj (j =t and I) and G|| for the inter-atomic and intra-atomic interac- 
tions, respectively. 

The Hamiltonian in Eq.Q can be further reduced by regarding all atoms as spin 
1/2 particles and defining collective angular momentum operators via Schwinger 
transformation, /,- = {btb_i + bi^bf)/2, Jy — {btbi^ — bi^bi)/2i, and J, = (btb^- — 

bibi)/2, alongside the raising and lowering operators 7+ = bfb:^ and J ~ b^bi. 
The reduced Hamiltonian is given by the expression, 

HDLMG=^fz+-7={a + a^)l+yj, (5) 



where we have dropped the constant energy term j[Ei^ + E^— 5 {G^ + G||)] + 
^ (G|| + G|| + 2G||). The average interaction energy of each atom is defined as 
^/N = j(G|| + G|| — 2G||) to account for the collective interaction among the 
condensates. The frequency detuning between cavity field and hyperfine transition 
is given by 4 = (»„ — coq. In the following analysis, only positive values of £, are 
considered (repulsive interaction) for the case of intra-species interaction, Gp^ or 
G^, larger than the inter-species interaction, G||. 

Equation (|5]l, the starting Hamiltonian for our analysis, corresponds to a general- 
ized Dicke Hamiltonian without the rotating wave approximation and an additional 
atom-atom quadratic nonlinear interaction. It is well known that the Dicke Hamilto- 
nian describes the collective dynamics for an ensemble of two-level systems driven 
by a quantum field cavity mode within a quantum electrodynamics (QED) configu- 
ration 1 14 1 . From many-body physics, the Hamiltonian in Eq. (5) is also equivalent 
to the Lipkin-Meshkov-GUck (LMG) model in the limit g ^0 U5J. Through the 
interaction among ensemble atoms, the LMG model, originally for A^ fermions dis- 
tributed in two A^-fold degenerate levels and interacting via a monopole-monopole 
force, was used to describe the Josephson effect in a two-species BEC and found to 
produce maximal pairwise entanglement of formation at the phase transition of its 
ground state [T6"l7l. Schemes proposed to implement a dissipative LMG model in 
optical cavity-QED 1 18 1 and in circuit-QED 1 19| have been discussed. The Hamil- 
tonian in Eq. (5) hereby will be called the Dicke-LMG (DLMG) Hamiltonian. Re- 
cently, it has been theorized that the ground state of the DLMG model supports 
phase transitions [20 1. In the semi-classical limit, both pitchfork and asymmetric 
bifurcations of the stable fixed points were found for this model II2TI . Also, the full 
quantum analysis of such a system showed a highly entangled ground state near the 
values of the semi-classical critical parameters ||22l . 

Now, we extend the DLMG Hamiltonian to include a Kerr medium ll23l l24l and 
degenerate parametric amplification |25| with nonlinear parameters K and j, re- 
spectively. The model Hamiltonian for such a system is given by the expression: 
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H = (Ooa^'d + K{a^af +x{a^+ a'^) + CO J. + |- J? + ^ (a + a') j\, (6) 

where the modified field frequency coq = (Of — K + 2% involves the modifications 

arising from the original Kerr and parametric amplifier terms, a^'^cP' and {a + a^) . 
Experimental realizations providing an assorted range of tunable parameters for 
the DLMG model may include a two-hyperfine-structure-defined-modes BEC cou- 
pled to a quantum cavity field mode through a one microwave photon process; e.g., 
trapped hyperfine ground states of a Sodium BEC inside a microwave cavity 1.261 . 
Arrays of interacting superconducting qubits coupled to the quantum field mode of 
a coplanar waveguide resonator may be considered as a physical realization limited 
by small ensemble sizes l.27il . 

In the literature, similar model Hamiltonians have been considered without the 
feedback from the BEC on the electromagnetic field; e.g., a second-order phase 
transition from immiscible to miscible is revealed in such a two-species systems by 
considering a linear mixing between the binary components |(28|. Stable domain- 
wall solutions, on top of flat continuous wave asymmetric bimodal states, can also 
be found near the point of the symmetry -breaking bifurcations |29J. The relative 
phase of domain-walls and breather-like dynamics for these dressed two-species 
BECs has also been studied by considering a classical external driving field If30ll31ll . 
Our results with a driven quantum field differs considerably from earlier studies: It 
is found that even in the semi -classical limit, the quantum field drive manifests in the 
excitation ratio parameter. When this excitation parameter is small or close to one, a 
discrete total excitation modifies the classical dynamics of the system strongly, pro- 
ducing localized asymmetric dynamics where some of the phase space trajectories 
present a running phase; e.g., Figph and Fig[3b. 

The proposed Hamiltonian in Eq. (6) may be experimentally realized by a two- 
species BEC, where the species are defined by two hyperfine-structure levels. The 
condensate is confined and coupled to a quantum cavity field mode through a one- 
microwave-photon process in the presence of a Kerr and a x medium. In general, 
the model described by Hamiltonian (J6]l has complex dynamics that deserve a study 
on their own. In order to give a clear interpretation, here, we limit ourselves to an 
analytical approach in the weak regime where intra-ensemble and ensemble-field 
couplings, as well as nonlinearities, are small compared to the field and hyperfine 
transition frequencies. 



3 Weak regime 

The full Hamiltonian in Eq. (|6]l does not conserve the total number of excitations, 
jV = 0^0 + J^, i.e., [H,.yV] 7^ 0. This fact prevents a simple approach like those 
given in the literature (|8] |2T] . This inconvenience may be bridged in the regime 
where nonlinearities and couplings are weak compared to the field and transition 
frequencies. In this weak regime it is possible to define a couple of unitary transfor- 
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mations ll32l 

f = el("'-"^'), (7) 

;/ =e-'K"+«'H. (8) 

The first of these transformations is equivalent to consider a squeezed basis; while 
the second one is similar to a polariton transformation and provides an effective 
rotating wave approximation, with the introduction of the small parameters Tj = 
X/oh < 1 and v ^ g{l- 2tj)/ (a)„ + a)o -4;^tj) <c 1, respectively. The following 
effective Hamiltonian, up to a constant and in units of h, is obtained by neglecting 
all the products of couplings with nonlinearities {i.e., K, ^, x^ ^^ ^ and v are all at 
least a couple orders of magnitude smaller than the field and atom frequencies) and 
moving into the frame defined by the total excitation number rotating at frequency 
(0 = 0^-4x7], 

Heff = 5f, + K (a'^a) ^+—jj + —— {aJ+ + a V"_ ) . (9) 



Here, the frequency detuning is given by 5 = (a^, — O and the effective ensemble- 
field coupling is A = (OV. This effective Hamiltonian in the weak regime is nothing 
else than the extended Dicke model studied in reference [22 \ plus a Kerr term. The 
Hamiltonian in Eq. (9) conserves the total number of excitations, [Hg//,^] = 0. 
Note that the exact dynamics defined by Hamiltonian Eq. (|9| can be calculated by 
quantum inverse methods 1.33 J . 



3.1 Semi-classical limit 

In order to fink our proposed configuration to a generalized pendulum problem, we 
apply the mean-field approach to study the semi-classical dynamics of this two- 
species BEC coupled to a optical cavity mode. It is possible to approximate the 
expectation values by considering the system in a separable state composed of a 
coherent photon state ||34l [35l . \^/ne'^), and a coherent spin 1/2 state 1361 l37ll . 
|z, 0), respectively. The expectation values for the field with the photon number n 
and optical phase are: 

(a) = a = ^/n e"^ , (10) 

(fl^> = a*=v«e-'^ (11) 

(fl^fl) = |a|2 = «. (12) 

If the ensemble is large enough, such that N^ I but less than the restriction brought 
by the two-mode approximation 1121 l38l . the expectation values for the ensemble 
operators may be approximated by those in the thermodynamic limit. 
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(/,) = -/z « f z, (13) 

(/±)^y±«^yr^e±'e, (14) 

where the fractional population difference is defined by the rotating angle in the 
corresponding atomic Bloch sphere, i.e., z = cosQ. The conserved quantity in our 
system is the mean total excitation number, which is given by the expression: 

(^)=^ = « + ^z. (15) 

This mean-field approximation plus the definition of a total phase variable, ^ — ^ + 
6, and an excitation ratio, k — 2yl^ /N, allow us to write the weak regime effective 
Hamiltonian in units of hNX/2 and up to a constant as: 

H= fA + jz]z+^2{k-z){l-z'^)cos<P, (16) 

where the re-scaled transition detuning now is shifted by the self-phase modulation 
from the Kerr nonlinearity, A — {5 + k)/X. Moreover, the re-scaled coupling ratio 
is defined as A = ^/X. This mean-field Hamiltonian is equivalent with that of the 
DLMG model in |22|, with the difference that in this case the Kerr and x non- 
linearities play important roles in the frequency detuning A and the characteristic 
interaction ratio A , respectively. 

Equation (16) may be viewed as a more general pendulum with nonlinear pendu- 
lum length and changeable pendulum mass described by the excitation ratio, k. This 
distinguish our system from the double-well configuration which is equivalent to a 
pendulum with just nonlinear pendulum length. Also, note the restriction k — z>0 
induced by the model Hamiltonian. Our system is equivalent to the case of a BEC in 
an asymmetric double well via a phase ;r-shift and a restriction given by ^ — z = 1 /2, 

Hadw^U + ^z) z-Vi^ coscp, (17) 



3.2 Fixed points of the system 

The dynamics of the semi-classical system is given by the equations of motions for 
the dimensionless fractional population difference and total phase variable. 



2{l-z^){k~z)smt>, (18) 

^ . . l+2fe + 3z2 

<P^A+Az COS0. (19) 

^2{l-z^){k-z) 
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Stationary states are found at the total phase variables values of <P ^0 and n. Due to 
the n phase difference with Eq. (fill, here the plasma and n oscillations will exchange 
places appearing at ^ = ;r and <P ~0, respectively. 

The fixed points of this Hamiltonian coincide with the critical points as z = 
dH/d<P and <P = dH/dz ETIl . Stationary states are found for the phase variable 
values <P —Q,7Z and the excitation parameter value 



k = 



3z2-l (I^z2)|(4+Az)| 



2z 



4^2 



A+Az\± iA+Azf~4z 



1/2 



(20) 



Notice that, in order to obtain a real excitation ratio, k, the allowed fractional 
population difference is bounded to the range z e [— l,z_] U [z+, 1], where z± = 
[2 - ZiA± 2(1 -ZiA)i/2]/A2 sets the condition 4 < 1/A. 

An example of fixed points in the weak regime is shown in Fig. l2] When the 
value of the excitation ratio k is less than or within the order of the magnitude of 
one, we have a quantum drive compared to a classical one (k ^ 1). It is the the 
introduction of the quantum drive and nonlinear processes that brings a peculiar 
breaking of the symmetry, different from the pitchfork bifurcation of the classical 
driving. It is possible to numerically sample the parameter space, {A,A,k}, and 
see that for any given excitation ratio, k, the fixed points satisfy mirror inversion at 
A —Q; i.e., z{A, \A\) — z(— A,— |A|). Also, if the frequency detuning, 5, is set to 
compensate the Kerr nonlinearity, 5 = — fc, it is possible to recover results that have 
been studied in the past; e.g. Fig.l2tb) andl2le) correspond to Fig. 2 in f2T|. 

From the mean-field Hamiltonian in Eq. ([T6|, it is straightforward to see that 
the quantum drive restricts the phase space accessible to Rabi oscillations for low 
excitation ratio. A: < 1, in order to keep the mean effective energy real, as shown 
in Fig. |3ja) and |3jb). Figures [3jc) and [3jd) allow us to see that the quantum drive 








10 



10 



10 

AAA 

Fig. 2 Example of fixed points in the weak regime, described by Eqs.i 18|l9|, 

(a-c) fc = 0.1 and (d-f) k = 10, respectively. Different tunning ratios, A = {S + k)/X, are shown: 

(a,d) A = —0.5, (b,e) A =0, (c,f) A = 0.5. Solid and dotted black lines correspond to stable and 



, for excitation ratios 



= 0.1 and(d-f)yt = 

= -0.5, (b,e) A = 

unstable fixed points for ;r— oscillations; while dashed grey lines correspond to fixed points for 
plasma oscillations. 
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also breaks the symmetry of symmetric initial conditions as isoenergetic lines are 
not symmetric with respect to the horizontal axis when the whole phase space is 
accessible. 

Figure HI shows an example of asymmetric Rabi and Josephson oscillations 
brought by the quantum driving field. A parameter set {A =0.5,k= 10} is taken and 
the fixed points found, Fig.|4|a). From the fixed points, two coupling ratios are cho- 
sen in the Rabi, A = 1, and Josephson, A = 8, regimes; this procedure delivers one 
stable fixed point, {A}, in the Rabi regime and two stable, B,C, and one unstable, 
{D}, fixed points in the Josephson regime. Figureplb) shows typical trajectories in 
these two regimes, where it is possible to see that symmetric initial conditions, i.e., 
{z±{t = 0) = ±ZQ,^{t = 0) = 0}, do not deliver symmetric trajectories. The latter 
can be seen straightforward from the position of the unstable fixed point, also called 
separatrix, and trajectories starting close to it. 

In the limit where the couplings ratio is large, A — > oo, it is possible to see a 
pitchfork bifurcations depending on the excitation ratio; an approximate critical ex- 
citation ratio can be calculated as kc w A^/2. At this critical value, a stable fixed 
point bifurcates into two new stable fixed points and the original fixed point be- 
comes an unstable fixed point acting as a separatrix in phase space. There is a sym- 
metry breaking in the dynamics due to the bifurcation, a transition from Rabi to 
Josephson oscillations; i.e., two initial symmetric states share identical dynamics in 
the Rabi regime, while in the Josephson regime they localize in different regions of 
phase space. Figure |5] shows an example of this symmetry breaking in the dynam- 
ics, in which for large coupling ratios, A ;:^ 1, i.e. 7] ^ A, it is possible to locate a 
pitchfork bifurcation point, kc+ w A^/2 even for the off -resonance condition, 8^0. 
Now the transition between the two-level system and the field 5 is different from 
zero but balances off the nonlinearity, i.e., A = 0. This condition, rj ^ X, relates to 



Fig. 3 (Color online) Nor- 
malized mean value of the ef- 
fective Hamiltonian energy in 
Eq. ^16^ , where a compression 
of the available phase space 
is shown for parameter sets 
{4,A,fe} = (a){-0.5,l,0.1} 
and(b) {0.5,8,0.1}. Also, 
the existence of asymmetric 
trajectories for symmetric 
initial parameters is intuited 
from the lack of symmetry 
with respect to the horizon- 
tal axis of the normalized 
mean value effective energy 
for parameter sets in the 
Rabi and Josephson regimes, 
{A,A,yt} = (c) {-0.5, 1,10} 
and (d) {0.5,8,10}, in that 
order. 
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a) 



10 




Fig. 4 (Color online) Fixed points and typical trajectories for parameter set {A = —0.5, A = 1 and 
8,^ = 10}. a) Fixed points for the equations of motion in the weak regime, described by Eqs. l |18| 
119), at ^ = showing the splitting of a single stable fixed point in the Rabi regime (solid purple) 
into two fixed points in the Josephson regime (dotted blue and dashed red). The unstable fixed 
point, or separatrix, in the Josephson regime (dash-dotted black) is also shown, b) Examples of 
trajectories with symmetrical initial conditions {z± = ±0.5, <P = 0}, leading to asymmetrical Rabi 
oscillations (purple dotted for z+ and dashed for z_) around the fixed point A; and asymmetrically 
localized Josephson oscillations around the fixed point B (blue dotted for z+) and C (red dashed for 
z+)- Two trajectories starting slightly above and below the separatrix D are also shown (dash-dotted 
black). 
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0, i.e. 5 
J18|19[, 


= -K 

at (P = 


a) Fixed points 
= in the large 



0.5 I -0.4 

k (units of A^) 

Fig. 5 (Color online) Off resonance, 5 7^ 0, with A = {S + k) /X = 
for the equations of motion in the weak regime, described by Eqs.l 

couplings ratio limit, A = i^/A — )■ 00 (numerical value A = 1000), showing the bifurcation of a 
single stable fixed point in the Rabi regime (solid purple) into two fixed points in the Josephson 
regime (dotted blue and dashed red). The unstable fixed point, or separatrix, in the Josephson 
regime (dash-dotted black) is also shown, b) Examples of trajectories with symmetrical initial 
conditions {z± = ±0.5, <P = 0}, identical Rabi oscillations (purple solid) around the fixed point A 
(numerical value k = 750000) and symmetrically localized Josephson oscillations around the fixed 
point B (blue dotted) and C (red dashed) Two trajectories starting slightly above and below the 
separatrix D are also shown (dash-dotted black) (numerical value k = 250000). 
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the phase space region where maximal shared bipartite concurrence in the atomic 
ensemble may be obtained in the quantum treatment of this model |22J. The differ- 
ence comes from the large excitation parameter ratio arising in this semi-classical 
analysis, kc+ » 1, i.e., Nq <S^n as Zc^ 3/A^ ^ 1. 



4 Conclusion 

In summary, we have presented an analysis of the classical dynamics of a two- 
species BEC large in size driven by a quantized field in the presence of nonlinear 
processes. In the weak regime, we find that the nonlinear phase from the Kerr nonlin- 
earity, K, shifts the transition detuning of the effective Hamiltonian of a generalized 
Dicke model; while the x nonlinear coefficient re-scales the coupling ratio. This 
mean-field Hamiltonian is equivalent to a nonrigid, nonlinear pendulum, for which 
a transition from Rabi to Josephson dynamics is identified depending on both the 
intra-BEC interactions to field-ensemble coupling ratio and the ratio between the 
total excitation number and the ensemble size. Moreover, we find that the symmetry 
of Josephson dynamics is broken by the quantum field, and an actual pitchfork bi- 
furcation point is found in the regime where the intra-ensemble interaction is larger 
than the field-ensemble coupling. Furthermore, It is known that symmetry-breaking 
in the classical dynamics of BEC may herald entangled quantum states |39 , 40|. Our 
results may provide a deeper understanding about the collective dynamics of inter- 
acting BECs and give another example in favour of the aforementioned conjecture 
relating classical and quantum regimes for nonlinear systems. 
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